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ABSTRACT 

The possibility of the magnetic flux expulsion from the Galaxy in the superbubble 
(SB) explosions, important for the a-U dynamo, is considered. Special emphasis is 
put on the investigation of the downsliding of the matter from the top of the shell 
formed by the SB explosion which is able to influence the kinematics of the shell. 
It is shown that either Galactic gravity or the development of the Rayleigh- Taylor 
instabilities in the shell, limit the SB expansion, thus, making impossible magnetic 
flux expulsion. The effect of the cosmic rays in the shell on the sliding is considered 
and it is shown that it is negligible compared to Galactic gravity. Thus, the question 
of possible mechanism of flux expulsion in the a-Q dynamo remains open. 
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1 INTRODUCTION 

The magnetic field in our Galaxy and in other spiral galaxies 
is usually believed to have been amplified from a weak seed 
field by a hydromagnetic dynamo, which exists due to the 
presence of the large-scale differential rotation and small- 
scale cyclonic turbulence in the Galaxy (Parker 1970, 1971; 
Vainshtein & Ruzmaikin 1971, 1972; Moffat 1978). ft has 
been suggested that any primordial magnetic field could be 
expelled from the Galaxy by the dynamic motions in less 
than a billion years (Parker 1971) so it would seem that 
some flux amplification is necessary to explain the Galactic 
field. The theory of such a dynamo has been formulated 
in a precise way through the mean field equations and the 
solutions to these equations indicate that the field would be 
amplified. The resulting field appears to correspond to the 
magnetic field patterns in our galaxy and others. 

On the other hand, a number of criticisms of this the- 
ory have emerged. One of them concerns the intense devel- 
opment of the small scale fields, which could damp the tur- 
bulence and stop the dynamo action unless they saturate at 
levels which do not interfere with mean field dynamo (Kul- 
srud & Anderson 1992; Parker 1992; Vainshtein & Cattaneo 
1992). 

Another unresolved problem is the expulsion of flux 
from the galactic disc. A very important point is that the 
theory of the a-Q dynamo predicts the amplification of some 
small preexisting magnetic field only if some magnetic dif- 
fusion is present in the Galaxy. However, too large a diffu- 
sion is destructive for the dynamo, because of its dissipative 
role in the process of generation and it seems attractive to 



suppose that a dynamo without diffusion at all will be the 
most effective. But Ruzmaikin, Shukurov, & Sokoloff (1988) 
showed that this is not possible. The physical reason for 
this is rooted in the very strong flux freezing of the galac- 
tic plasma, because magnetic lines cannot break and the 
number of field lines in the disc can be increased only by 
toroidal stretching, which is accomplished by the dynamo 
action. Any stretching creates field of both signs, to con- 
serve the total flux, and, for net amplification to occur, those 
portions of field lines which are of the wrong sign must be 
expelled from the Galaxy. In the standard a-Q dynamo the- 
ory this is done by magnetic diffusion. 

Ruzmaikin, Shukurov, & Sokoloff (1988) considered 
a-f2-dynamo in the case of the thin disc and demonstrated 
that temporal evolution of the magnetic flux is governed by 
the following set of equations: 



+ G / B r dz, 



(1) 



(2) 



where G = rdQ/dr-measure of the differential rotation, /3- 
magnetic diffusion, and and 1 correspond to the center 
plane and boundary of the Galactic disc. 

The j3 terms represent the expulsion of flux. It is clear, 
that setting /3 = in (UJ) and (0) we immediately get that 
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Brdz = const, 



(3) 



and that J B^dz may grow only linearly due to the stretch- 
o 

ing of lines by the galactic differential rotation. Thus, there 
is no exponential field growth, essential to the dynamo. One 
can easily see that for the dynamo to operate, there must 
be some nonzero flux escape through the upper boundary of 
the disc, that is /3 0. 

This is consistent with the topological constraint that 
the total number of lines of force including those negative 
lines expelled from the disc must be constant, fndeed, a sim- 
ple estimate of the expulsion terms making use of the numer- 
ical results in Ruzmaikin, Shukurov, & Sokoloff (f 988) shows 
that the negative flux expelled during one growth cycle is 
comparable to the positive flux in the disc at the beginning 
of the e-folding. 

Although the a-Sl dynamo theory is complicated, this 
physical intuition of flux expulsion can be considered in the 
absence of these complications. Further, the mechanism of 
expulsion need not be tied to the (3 diffusion inside the disc. 

The main problem with the expulsion of flux is that 
this flux is loaded with matter so that it is related to the 
expulsion of matter against the strong gravity of the galactic 
disc. The most likely process to expel flux is the phenomenon 
of sequential supernova (SN) or superbubbles (SB), which 
sweep up matter into dense, radiatively cooled shells. Mag- 
netic field, tied to the matter due to the strong flux freezing 
in the ISM, is also swept up and deposited in these shells. 
If some part of the shell leaves the Galaxy, it carries the 
frozen- in magnetic field with it, thus producing the flux ex- 
pulsion. 

But most of these superbubbles are not powerful enough 
to expel matter out of the gravitational well of the disc. The 
only possibility for flux expulsion seems to be: as the bubble 
expands, the field lines in the shells of SBs are not horizontal 
but form arcs, along which matter can slide down, lowering 
the amount of matter on the top of the lines and allowing 
some flux to escape. 

There is further difficulty with the mechanism involv- 
ing the SBs which is relevant to its application to the a-Cl- 
dynamo. The a-f2-dynamo assumes small-scale turbulence 
while the cavities produced by SB explosions may reach 
~ 500 pc or larger, which is greater than some of the length 
scales of the galactic disk. For the dynamo theory in its con- 
ventional form to be applicable it is important that turbu- 
lence be small scale, because it involves the expansion of the 
turbulent electromotive force £, which describes the effect 
of turbulent motions on the mean (or ensemble- averaged) 
magnetic field, in terms of the mean magnetic field itself and 
its spatial derivatives: 



£i — a-ij < Bj > +f3, 



d < Bj > 



ijk " 



d. 



(4) 



(Moffat 1978) 

If the scale of the turbulence is too large, then the ex- 
pansion is invalid and usual a-il-dynamo theory must be 
modified. That's why, for example, direct application of a 
and /3 tensors calculated by Ferriere (1995, 1998) for SBs 



and SNs to the a-fi-dynamo theory can lead to an overly 
optimistic estimates of the rate of flux escape. 

Although a-fi-dynamo theory is not strictly applicable 
to the case of SBs, the actual operation of them in amplifying 
the field is clear from the work of Ferriere (1991, 1995, 1998). 
It is also clear that the rapid escape of the flux from the disk 
is essential. In this paper we show that because of the deep 
gravitational well of the Galaxy it is difficult for the matter 
and field lines to escape, and consequently for the mean field 
to grow. In Ferriere's works she finds the lines of force rising 
with the SB but does not follow them long enough to see 
that they must fall back into the disk and inhibit the growth 
of the field. 

In this note we quantitatively examine the dynamics of 
the rising field lines carried by SBs and show that even with 
sliding the matter and flux are unlikely to escape. Thus, 
the requirement of the escape of flux provides a strong con- 
straint for the a-Q dynamo to overcome if it is to amplify 
the galactic magnetic field. 



2 SLIDING OF MATTER FROM THE TOP OF 
SB; FORMULATION OF THE PROBLEM 

Multiple supernovae from OB associations can carve out 
large cavities of hot gas, called superbubbles. McCray & 
Snow (1979) first described them. When the SB expands 
into surrounding medium it sweeps up interstellar matter, 
giving rise to a massive expanding shell. Inside the volume 
surrounded by this shell a hot rarefied low-density gas is con- 
tained which provides the pressure driving further expansion 
of the shell. The energetic source for sustaining this pressure 
is provided by continuous energy input from the SN explo- 
sions in the center of SB. Weaver et al. (1977) calculated 
the evolution of the bubble driven by the continuous wind 
from the central source and Mac Low & McCray (1988) ap- 
plied this theory to the case of supernovae driven SBs. They 
show that the radius of such a SB, expanding in a uniform 
medium of density po with continuous energy input in the 
center Lsn = LzgW 38 ergs s _1 (the luminosity of SB con- 
veniently expressed in units of I0 38 ergs s _1 ), is given by 



R(t) 



( 125 y /5 L^ - i/5 t^ 



267 



-^38^7 

n 



1/5 



pc, (5) 



where tj = t/l 7 yr, with the velocity of the envelope chang- 
ing as 

1/5 



L38 

n t? 



u(t) = R(t) « 15.7 [ ) km s _1 . (6) 



The inner pressure in the volume bounded by the shell varies 
as 

P 01- 7 r 2/5 3/5,-4/5 

nn[t >- (3850tt)2/5 SivPo 



4.1 x 10 



r2 3 \ V 5 
-12 / ^38 n \ -3 

ergs sm , 



'-7 



(7) 



due to the work done on the expansion of the shell and the 
energy injection in the center of SB. 

When the shell expands in the real Galactic environ- 
ment, there is also a gravitational force which tends to 
slow down the vertical expansion. Also, the distribution of 
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medium is highly inhomogeneous with height z over the 
Galactic plane. Various components of ISM have different 
length scales and characteristic densities, but they all have 
exponential or Gaussian decreasing profiles in z and thus 
drop very rapidly with height. This effect is very noticeable 
for powerful superbubbles, for which the expansion radius 
may exceed the height scale of the matter distribution. As 
we will see this can lead to Ray leigh- Taylor instabilities. 

At the final stage of the SB expansion, its velocity be- 
comes comparable to the sound velocity of the ambient ISM 
and shock wave will no longer exist. However, this does not 
prevent the SB from further expansion, because there is still 
a lot of momentum in its massive shell of swept up material 
and the ram pressure at high z is negligible. This means that 
expansion actually continues until the velocity of the shell 
drops to zero under the action of Galactic gravity: 



u = 0. 



(8) 



This stopping condition is different from the one which 
Ferriere used (1998) in her investigation of the role of SBs 
in the a — Q. dynamo. Her condition (expansion stops when 
the shell velocity is of the order of speed of sound in the 
ambient ISM, in other words when the pressure inside the 
SB becomes comparable to the ram pressure) is applicable 
only for massless shells without any inertia. In reality shell is 
quite massive and the momentum stored in it drives further 
expansion of the bubble against the ram pressure and the 
galactic gravity. 

Also when the SB expands in the presence of the gravity, 
the swept-up matter in its shell is likely to slide downward 
along the shell as mentioned before. This can influence the 
efficiency of flux removal by SBs. Obviously, in the absence 
of sliding, escaping matter takes with it all the frozen-in 
magnetic flux. Now, if we allow the matter in the shell to 
slide downward perpendicular to the field B in the shell it 
will take the magnetic field lines away from the top of the 
SB, thus reducing the magnetic flux to be removed. But 
it can also have a significant effect on the dynamics of the 
SB itself, because as the matter slides from the top in any 
direction, the upper part of the shell becomes lighter and 
acceleration due to the inner pressure of the hot gas will gets 
larger in proportion to the surface density decrease, while 
the gravitational deceleration stays the same. This results in 
some additional acceleration of the top of the shell and, in 
principle, this can significantly change the dynamics of this 
part of SB if there is enough time before the shell stops. If 
this happens, and a substantial part of the mass slides down, 
the top of the shell may continue its expansion upward and 
drive the remaining matter and frozen-in flux to a further 
distance from the galactic plane and possibly expel the flux 
entirely from the Galaxy. For this reason it is very important 
to estimate the numerical value of this effect. 



tion is restricted to the case when the angular distance of 
the shell element from the shell top, 9, is very small, 9< 1. 

We suppose the expansion of the SB is described by 
some law R = R(t). We will also suppose for simplicity that 
the form of the shell near its top can be roughly approxi- 
mated by a sphere. The validity of this assumption will be 
discussed later. 

We made some further simplifications, one of which is 
the neglect of the pressure gradient along the shell. It is clear 
enough that any such pressure gradient would only reduce 
the downsliding of the matter. This means that our estimate 
is only an upper limit of the sliding and the real sliding will 
actually be smaller. 

We also completely neglect the influence of the mag- 
netic field on the dynamics of the sliding. This assumption 
is justified in the early periods of the Galaxy's life, when 
the magnetic field was weak. At the present time this is 
not completely valid, because the magnetic field is strong 
and magnetic tension may play some dynamical role in the 
expansion process. SBs with strong magnetic field were con- 
sidered analytically by Ferriere (1991) and numerically by 
Tomisaka (1992). 

Thus, the only external force in our analysis is the grav- 
ity due to the stars and ISM in the Galaxy. We take gravity 
as given, because the self-gravity of the bubble is negligible. 
With this in mind the system of the equations of Giuliani 
describing the gravitational fall of matter from the top of 
the shell reduce to 



^(R 2 sin 9a) 



pouR 2 sin 9 — (R(t) sin 9ai 

Uu 



or 



+ 



V II 

+ i 



dvu 



0. 



(9) 



(10) 



Here i>|| is the tangential velocity of the matter along 
the shell arising from the presence of gravity, u — R(t) is 
the velocity of the shell expansion (dot means time deriva- 
tive), a is the surface density of the shell, po is the density 
of the unperturbed gas in front of the shell, and g is the 
gravitational acceleration. 

Now, at an early stage of SB expansion, when its ra- 
dial velocity is very high and the transverse velocity due to 
the galactic gravity is not very large, the effect of sliding 
is negligible, because the total velocity of the shell element 
is directed almost radially. Thus, it is reasonable to suppose 
that the most noticeable effect of sliding will occur only dur- 
ing the later stages of the shell evolution, when it sufficiently 
slows down in the radial direction. But by the onset of this 
later stage the SB has expanded in the direction perpen- 
dicular to the galactic plane to a distance larger than the 
scale height of matter distribution. This reasoning allows us 
to neglect the second term proportional the ambient density 
in equation ( |lo| ) . It makes equation ( |lo| ) completely inde- 
pendent of equation (0) since it then contains no terms in 



3 BASIC EQUATIONS 

Giuliani (1982) gave a general formulation of the thin-shell 
approximation for hypersonic, hydromagnetic flows, axisym- 
metric about the z axis, including motions along the shell. 
We use his equations in our description of gravitational mat- 
ter removal from the top of SB. All our further considera- 



However, in equation (|9|), we must keep terms in po 
because a still grows due to swept up matter. 

Also we can neglect the term dv\\/d9 in ( |lo[ ) compared 
to the radial velocity u. We discuss this omission later. Thus 
equations (|^) and (|l^) reduce to the following simplified 
system of equations: 

(R 2 (t) sin 9a) = p uR 2 (t) sin 9 
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(11) 



(12) 



dv« u . „ 

In general, we assume that gravitational acceleration, 
g, is a function of the z-coordinate in the Galaxy. Bearing 
in mind that u = R(t) and v\\ = at t = 0, we can integrate 
equation dl2|) to get 



sin 



R(t')g(z + R(t'))dt', 



(13) 



where zq is the height in the Galactic plane where the ex- 
plosion occured. 

We see that vn oc sin 9. This means that in equation dll| 



we may neglect the last term near the top of SB since it is 
proportional to sin 2 8. Then equation ( |l l| ) further reduces 



to 



R 



2R + 



d 



sin 9 89 



sin 0v\\ 



poR- 



(14) 



This equation can be integrated with the initial condition 
a = at t = 0: 



*{t) 



R 2 (t) 



exp 



K(t')dt' 



x J R 2 Rp(z + R(t'))exp I J K(t")dt" J dt', (15) 



o 

where 

sm9 09 v 

Substitution of expression (113) into (hq) gives 



(16) 



K = 2 cos f 



1 



R(t')g(t')dt' 



R 2 (t) 



R{t')g(z + R(t'))dt', 



(17) 



near the top of SB. 

From the formulae ( |l5| ) and ( p^ ) it is easy to see that 
the importance of sliding near the top is determined by the 
quantity 



o 



cm 



With slidin 




o 



PS, 



Without sliding 



N = 75 



t, Myr 



With slidin 




Without sliding 



40 60 

Number of SNs 

Figure 1. (a) Dependence of polar radius R (at = 0) upon the 
expansion time t for the case of the sliding of matter due to the 
Galactic gravity (upper curve) and without sliding (lower curve). 
The case of 75 SNs in SB is considered, corresponding to the 
luminosity L = 7.5 X 10 37 ergs s — . Note that SB with sliding has 
larger final size, than that without sliding, (b) Dependence of the 
maximum radii of the SB, Rmax, on the luminosity (exp ressed 
through the number of SNs in the SB in equation fell)) with 
(upper curve) and without sliding(lower curve). The difference in 
final sizes of the SBs with and without this effect increases with 
the bubble's luminosity L. 



= e 



exp 



R 2 (t) 



R(t')g(t')dt' 



(18) 



which can be determined for any given R(t). 

If £ ~ 1, we can safely neglect the sliding of matter but 
if £ 2> 1, sliding will play important role in the dynamics of 
late stages of SB expansion. 



4 INFLUENCE OF SLIDING ON THE SHELL 
EXPANSION 

Let us examine the equation for R(t) near the top of the 
bubble taking into account the effect of the unloading the 
matter from the shell's top on the radial expansion of the 
SB itself. 

To do this we consider a solid angle dfi of the shell near 
its top. We can write the following equation for its motion: 



© 2000 RAS, MNRAS 000, [|-?? 



Magnetic blux Expulsion in the rowerjul buperbubble Explosions and the a-il Dynamo. 5 



Tit 



(aRR 2 dQ.) = (P in - P out ) R 2 dQ - R 2 agdQ. 



(19) 



We can combine this equation with equation ( J15[ ) to find 
R(t) and a(t) near 6 — as a functions of time t if we know 
how Pi n behaves. All other quantities, namely Pout and g 
are given empirically, as a function of z, by the position of 
the shell's top. 

Equations ( |l9| ) and ( Jl4| ) can be combined to give 



R + R 



2P0 



2R 
~R? 



R{t')g(t')dt' 



p. 

1 in 



(20) 



and together with equation (^) give R(t). 

These equations are solved numerically to get the be- 
havior of R. First we supposed that the inner pressure is 
governed by equation (Q). This underestimates the pressure 
and the height reached. But if there is no escape of flux 
under this assumption there is certainly no escape in real 
conditions. We have taken the model of ISM from Ferriere 
(1998), that is we supposed that density and pressure of 
ISM are contributed by 5 components, having different num- 
ber densities, length scales and temperatures: neutral, cold, 
warm, ionized, and hot. We also use her approximation for 
the gravitational acceleration g. We suppose for simplicity 
that the luminosity of the SB is constant in time during 37 
Myr, until the death of the 8Mq stars (Ferriere 1995), and 
equal to 

L = 10 36 x N erg s"\ (21) 

where TV is the number of SNs in the star cluster. For this 
calculation we suppose that the inner pressure Pi n changes 
in accordance with equation (f?|) during the first 37 Myr of 
the SB expansion, with po the mass density at the site of 
explosion. After 37 Myr the interior cools adiabatically be- 
cause the inner pressure and inertia of the shell continue to 
drive the shell expansion so that the volume bounded by the 
shell increases. We also carry out calculations for the differ- 
ent pressure law. But we neglect radiative cooling of the hot 
gas in the interior of the bubble during the entire explosion, 
so that real SB expansion is always smaller than we obtain 
here. 

In Figure [l|a we show the dependence of R upon time t 
for the case of a SB with N — 75 SNs in it, going off near the 
Sun at an initial Galactic altitude 100 pc. For comparison we 
also depict the curve without sliding for the same SB, which 
is obtained by setting g in the integral in the left-hand side 
of equation (|20|) to zero (but not in the right-hand side!). 

If the sliding is taken into account, the time when ex- 
pansion stops is t s = 35 Myr and Rmax = 1336 pc; neglect- 
ing sliding we get t s = 30 Myr and R ma x = 1094 pc. It 
is obvious that effect of sliding should be more pronounced 
in more powerful bubbles because of their longer lifetimes 
and stronger gravity at the heights to which they bring the 
matter. To illustrate this we plot in Figure Qb the maximum 
radii of SBs of various luminosities with and without sliding 
for the same conditions as the Figure [I]a. It is clearly seen 
that sliding plays important role for powerful SBs, making 
their final radii dozens of percent larger than that without 
sliding. 

In Figure ^| we plot, for comparison, the expansion ve- 
locity of the N = 75 SB top and the velocity of sliding along 
the shell for the same SB, divided by sin 8 (which virtually 



>" 




Figure 2. Plots of the expansion velocity of the SB top u (de- 
clining curve) and the velocity of the downsliding of the matter 
along the shell (growing curve), divided by sin0. These plots 
show that at the early stages of the shell expansion we can neglect 
dv\\/89 compared to u in equation 



equals to dv\\/d9 near the SB top). One can see that, up 
to the first 20 Myr, u is larger than dv\\/d8, which justifies 
our neglect of corresponding term in equation ([tC)|). By that 
time most of the SB expansion has already occured, so the 
inclusion of the term with dv\\/d9 into ( fl2|) does not change 
the final results significantly. Moreover, if we do include it, 
it would only suppress sliding, as can be seen from @, 
so that our results for sliding and shell expansion can be 
considered to give upper bounds for the final height of the 
shell. 



5 RAYLEIGH- TAYLOR INSTABILITY 

It is tempting to expect that for even more powerful SBs, 
than JV = 100 SNs, the effect of sliding will be even more 
pronounced and the final size of the SB may reach scales 
comparable to the size of Galaxy, thus making possible the 
escape of the flux. But when the luminosity of the SB ap- 
proaches 10 38 ergs s _1 , corresponding to the number of SNs 
TV ~ 100, another effect becomes important for the fate of 
SB. At such a large luminosity the shell starts accelerat- 
ing some time before its expansion could be stopped by the 
gravity and it would accelerate, in principle, to a very high 
velocity if there is enough time for it. 

However, it was first noticed by Mac Low & McCray 
(1988) and then proven numerically by them (Mac Low & 
McCray 1989) that as soon as the shell starts to accelerate 
it becomes Rayleigh- Taylor unstable and eventually breaks 
up. They called this process a "blowout" of the shell into 
the halo. Indeed, the effective gravity in the moving shell in 
the case of acceleration is directed towards the center of SB, 
that is the dense cold gas in the shell is pushed by rarefied 
hot gas of the interior and this leads to the instability. The 
shell fragments into blobs of cold, dense gas, which continue 
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moving with velocities they attained before fragmentation. 
There is no further significant acceleration of the shell, be- 
cause the hot rarefied gas from the interior of the SB escapes 
into halo thus allowing the inner pressure to drop. After 
that time each blob will move ballistically, sweeping up some 
small mass, though this effect might be not so important at 
high altitudes, where the density of all the components of 
the ISM is very low. This means that if the speed of the blob 
at this time is less than the escape velocity from the Galaxy, 
matter can not leave but must return to the Galactic plane 
from halo and there is no contribution to the flux escape. 
For this instability the galactic gravity contributes to the ef- 
fective gravity as well, thus making the shell unstable even 
when it is still decelerating. 

In our situation, when we have a semi-infinite hot rar- 
efied medium and a dense slab of gas of finite thickness H up 
on it, the fastest growing mode of the instability has a scale 
of the order of H , so that the growth rate of the instability 
is given by 



7 2 = - (R + g) /H. 



(22) 



The increment of instability or the amount of e-foldings re- 
duces to 



~/dt = / dt^J (g + R) /H. (23) 
For the fragmentation to proceed effectively we require 



jdt > 1. 



(24) 



To get the shell thickness, let us note that the gas ini- 
tially located between the heights zo + z and zq + z + dz is 
deposited into the shell between h + dh and h from the shell 
outer surface. The conservation of the number of particles 
accounting for the sliding gives that number density of par- 
ticles in the shell at local thickness h related to the number 
density of particles initially at the point zo + z as 



— no(zo + z)z dz 

(To 

where 



n{h)R dh, 



<T0(t) = 



R?{t) 



R 2 Rp(t')dt', 



(25) 



(26) 



is just a without sliding. The number density n is given 
by n = P s h(h)/kT s h, where T s h is the temperature in the 
shell, which we assume to be equal to 10 4 K, and P B h is 
the pressure at the local thickness h in the shell. Here we 
have taken into account only the thermal pressure and ne- 
glected the cosmic ray pressure. This seems to be the rea- 
sonable assumption, because, as will be shown later, they 
have very high drift velocity and may easily escape from the 
compressed shell along the magnetic field lines deposited 
into it, since the lines themselves leave the shell. The pres- 
sure P s h is comparable to the inner pressure Pi n , because 
inside the shell it has to drop from P;„ on the inner surface 
to Pout <C Pin outside the shell. 

Combining all these considerations we obtain that 



(To , 



n (z + z)kT sh z' 
P n R 2 



dz. 



(27) 



o 

Dh 



Stopped SBs 




RT-unstable SBs 



200 400 600 800 

Number of SNs 



Stopped SBs 




RT-unstable SBs 



Number of SNs 

Figure 3. (a) Final size of the SB, R m ax, as a function of N 
- the number of SNs producing different luminosities. Shells of 
the low luminosity bubbles are stopped by galactic gravity (left 
branch), while those of the more powerful bubbles (right branch) 
are disrupted by the Rayleigh- Taylor instability, (b) Distribution 
of times when the SB expansion stops, t max , either due to the 
galactic gravity (left branch) or the Rayleigh- Taylor instability 
(right branch) for different number of SNs. 



Numerical estimates show that during the first several 
e-foldings of the Rayleigh- Taylor instability, the shell thick- 
ness H changes very little, less than 10%. During this time 
the size of the shell changes less than ~ 30%, so that the 
geometrical effect of the stretching the scale of the pertur- 
bation mode in the expanding shell is quite moderate. For 
this reason the modes which were unstable at the very on- 
set of the instability stay unstable during several e-foldings 
thus developing the nonlinear stage of the instability and 
disrupting the shell. 

During the development of the instability the velocity of 
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the shell does not change drastically, it is larger only (3—5)% 
than the velocity at the very onset of instability. Thus, we 
may safely assume that after the Rayleigh- Taylor instability 
is fully developed in the shell, we get no further acceleration 
of the shell fragments. 

We consider the role of the Rayleigh- Taylor instability 
on the fate of SBs of various luminosities with inner pres- 
sure from (Q) located at the altitude zq = 100 pc near the 
Sun. The results for maximum size of the shell and time 
when it either stops due to gravity or fragments because of 
Rayleigh- Taylor instability, are shown on Figure |S| Equation 
(M) is chosen to be the condition for Rayleigh- Taylor frag- 
mentation of the shell, because it corresponds to the onset 
of the nonlinear stage of this instability when the shell is be- 
ing disrupted. The two branches of Figure ^ correspond to 
the bubbles which were stopped (left branch) and to those 
which were disrupted (right branch). We see that for a cho- 
sen dependence of inner pressure upon time the transition to 
the Rayleigh- Taylor regime occurs for N = 104 supernovae 
in the SB. The maximum possible size and the greatest life- 
time of the SB are achieved just before this N and are equal 
to Rmax = 2616 pc and t max — 52.6 Myr. After that size is 
reached the expansion drops rapidly due to the early onset 
of the Rayleigh- Taylor instability in the shell, and we see 
that for SBs with N ~ 1000 shell fragmentation occurs at 
a very early time ~ 5 — 6 Myr. It must be emphasized that 
at such an early stages the approximation of constant lu- 
minosity may not be justified and dynamics might be more 
complex. We believe that the result of the analysis is not 
very sensitive to the model assumptions. 

In Figure ^ we plot the dependence of the blob veloc- 
ity upon the luminosity of the SB. One can see that it is 
sufficiently less than ~ 430 km s , the estimated lower 
bound on the Galactic escape velocity (Leonard & Tremaine 
1990), so that we may conclude that shells formed by SBs 
with number of SNs in them N < 10 3 do not give rise to 
the mass and flux outflow from the Galaxy. Even if we go 
to a SB with a luminosity an order of magnitude larger 
(N — 10 4 ), the velocity at the moment of fragmentation 
is only Vmax ~ 336 km s _1 , which is obviously not enough 
to leave the Galaxy. 



6 IMPORTANCE OF THE SHAPE OF THE SB 
TOP 

In our treatment of the SB expansion we have considered 
the shape of the shell near its top to be spherical. This en- 
ables us to use a simple fact that in this case the projection 
of Galactic gravitational acceleration along the shell is just 
gu — g sin 6 which sufficiently simplifies the problem. 

In reality, of course, the surface density is not uniform 
on the top but depends upon the angle 8. The nonunifor- 
mity grows with growing 9. The inner pressure will acceler- 
ate parts of the shell closer to the top stronger than ones 
further from the top and it will distort the form of the shell. 
This distortion in its turn changes gu which influences the 
sliding of the matter and thus leads to further changes of the 
shell's form. The accurate treatment of the problem requires 
including this effect self-consistently in our calculations, but 
we can avoid this by noting that influence of a change of the 
shell's shape on the process of sliding can be attributed to 




400 600 

Number of SNs 

Figure 4. Velocity of blobs formed in the process of the frag- 
mentation of SB shell, v m ax, versus the luminosity of the SB 
(number of SNs in the initial star cluster). Though this velocity 
increases with the luminosity of the SB, it does not reach the 
escape velocity from the Galaxy. 



the change in gravitational acceleration g, rather than the 
projection angle. 

We carried out calculations identical to those with 
spherical top but have taken g in equation (^) to be 4 times 
larger than it is in reality. The result was that, as the lumi- 
nosity of the SB was increased, they expanded faster, due to 
the more effective sliding of the matter from the top. But this 
in turn lead to a more rapid onset of Rayleigh- Taylor insta- 
bility in the shell: it started to develop for SBs with N > 69. 
The bubble with iV = 69 stops at a time tar = 45.7 Myr and 
size Rmax = 2398 pc. The conclusion is obvious: the change 
of the shell shape might influence the sliding of the matter, 
but it leads to the development of the Rayleigh- Taylor in- 
stability for even smaller luminosities than in the case of the 
bubbles with the spherical top, making the impossibility of 
the expulsion of the flux from Galactic disc even more cer- 
tain. That is why we think a more self-consistent approach 
to the problem of the shell shape will not change the general 
result. 



7 DIFFERENT PRESSURE LAW 

The pressure law (0) which we used in all our calculations 
was derived actually for the case of the SB expanding in 
a uniform medium and thus may be not a very good ap- 
proximation for our purposes especially when the sliding of 
the matter influences the expansion of the shell and its size 
cannot be described by equation ^ . For that reason we de- 
cided to use a different, more realistic pressure dependence 
to check if it makes a significant difference in our results. 

Maciejewski & Cox (1999) proposed a simple, explicit, 
analytical approximation for the kinematics of the blast 
wave propagating in an exponentially stratified medium: 
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po = P*e 



z//i 



(28) 



where p* is the density at the explosion site and h is the 
stratification length scale. They considered an explosion in 
the framework of the Kompaneets approximation (Kompa- 
neets 1960) with inner pressure constant throughout the vol- 
ume engulfed by the shock. They showed that the form of 
the shell is very close to an ellipsoid with minor and major 
axes b and a related by 

tan -4- = sinh (29) 
2h 2h v ' 

and the distance from the explosion site to the ellipsoid cen- 
ter, s, 

S d 

tan — = cosh — . (30) 
2h 2h v ' 

At the same time, Weaver et al (1977) showed that in- 
ternal energy of the SB interior is constant fraction of the 
total energy and for a spherical SB expanding in a homoge- 
neous medium 



o 
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(31) 



We assume that this is also valid for the case of SB in a non- 
homogeneous density distribution, so that that inner pres- 
sure 

Ei n _ 5 LsNt 

~V~ ~ IT nab 2 ' 



P — 

i Ml 



(32) 



We relate the semi-major axis a to the distance from the 
center of explosion to the top of the shell: a + s = R, by 



a ( a \ R 

h log cosh — — — 

2h 6 I 2hJ 2h 



(33) 



Formulae ([B2j),([33j), and ( [29] ) give us P in for a given R. 
The approach of Maciejewski & Cox (1999) includes neither 
galactic gravity nor the slippage from the top but it takes 
the inhomogeneity of the surrounding medium into account 
and enables us to test the stability of our results against 
different model assumptions. 

The real ISM contains many components distributed 
with various length scales so we take rather arbitrarily h = 
200 pc in our case. The particular choice of h turns out not 
to play a significant role. The results for superbubbles of 
various luminosities located at zq — 200 pc near the Sun are 
shown in Figures ^a and ^b. 

We see that differences are only quantitative compared 
to the case of pressure law ([?]). Strong Ray leigh- Taylor insta- 
bility starts to dominate the kinematics of the shell when the 
number of SNs in SB is larger than N = 18. This SB reaches 
the maximum size R ma x = 447 pc at a time t max = 41 Myr. 
Due to the specifics of the the chosen pressure law, SB with 
the higher luminosity, developing the Rayleigh- Taylor insta- 
bility, reach somewhat larger size, maximum is Rrt = 630 
pc. If we take, for example, h — 100 pc, then the shells 
start to be not stopped by the gravity but disrupted by 
Rayleigh- Taylor instability even for smaller N. Thus, again, 
for powerful SBs, expansion is limited by the shell fragmen- 
tation so that all major results of the consideration with the 
simplified pressure law remain valid. 




400 600 

Number of SNs 



Figure 5. (a) Final size of the shell, R m ax, and (b) velocity 
of the shell, v max , at the moment of stopping or fragmentation 
versus the number of SNs in SB. Pressure law is different from 
the case represented on the Figure ^ and is given by equations 
and (§§). 



8 POSSIBLE IMPORTANCE OF THE CR 
PRESSURE IN THE SHELL 

Kulsrud (1999) proposed that sliding of the matter from 
the top of the SB shell may be inhibited to some extent by 
cosmic ray (CR) pressure gradient, thus further supporting 
our conclusion about the impossibility of flux expulsion from 
the Galaxy. Now, on the basis of the better knowledge of the 
processes going on in SB shell, we can check this idea in more 
detail. 

Let us consider the magnetic flux tube with the cross 
section constant along the tube in the shell which reached 
the size R and expands with velocity u — R. This assump- 
tion is good for the cylindrical shell, with its axis lying in 
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the Galactic plane, but it will be clear further that this as- 
sumption is not important. The continuity equation for the 
CR along the flux tube reads: 



d_ 



(n CR R) + -kz [ 



ncR [Vi 



+ v d )] 



0. 



(34) 



where ucr is the number density of the CR, vu and Vd are 
the velocities of matter along the shell with respect to the 
rest system and of CR with respect to the matter corre- 
spondingly. The drift of the CR along the magnetic field 
is determined by the gradient of their number density and 
by the scattering of the CR by the Alfvenic turbulence. The 
scattering of CR by the self-generated Alfven waves was first 
considered by Kulsrud & Pearce (1969) and Wentzel (1969) 
and they showed that CR cause the Alfven waves to grow 
with a rate 



r 



„tt v d - Va n C R{e) 

O— Ho , 

4 va n 



(35) 



where Qo is the nonrelativistic cyclotron frequency, nci?(e) 
is the number density of CR with energy greater than e, n 
is the density of the ISM, va = B /^/innrriH is the Alfven 
speed (niH is the mass of hydrogen atom), and C is a con- 
stant of the order of unity whose value depends upon the 
energy spectrum of CR. These waves are also nonlinearly 
damped by the scattering of the beat waves on the ions: 

(Lee & Volk 1973; Kulsrud 1978) with SB being the mag- 
netic field perturbations due to the Alfven waves, Q - rela- 
tivistic Larmor frequency, and vt - thermal velocity of ions 
in ISM. Scattering of CR occurs at a rate 

- = «(f) 2 , (37) 

(Kulsrud 1995), so that the typical mean free path of the 
CR is 

On the other hand, the drift velocity is given by 

A dn C R , . 

ncRVd = c- — . (39) 

Assuming that Vd ^ va, that Alfven waves are in equilib- 
rium, r = jd, and combining equations (^),(|^), (^), and 
© we get 

1/2 

' 2 cvtUVa 



"d 




= A 



1 dn C R 
R 89 

dnci 



Rn 2 CR 



(10 



1/2 



(40) 



If we introduce the characteristic drift velocity 



ucrR 
and suppose that 

VD > W||, 



(42) 

then it is easy to see from equation (^) that when Vd ~ v\\ 
the variation of the CR density can be expressed as 
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Figure 6. Characteristic drift velocity of CR, vu, in the shell of 
SB due to the scattering of the CR by the self-generated Alfven 
turbulence. Note the significant difference in the magnitudes of 
vd and the sliding velocity vu in the shell (Figure H). 



(43) 



VD 



ucrR = n C R Q Ro ^1 + (^jj) Xi(#)^ . 

where Xi($) ' s some function of order of unity. 

To see if the condition ([I2]) is fulfilled, we calculated ' 
given by ([ll]) and (^) at various time moments for the SB 
with TV = 75 SNs, taking into account the compression of the 
matter and CR in the shell, which we calculate in the manner 
similar to the calculation of the thickness of the shell, and 
supposing that the temperature in the shell after cooling 
is T s = 10 4 K. Compression of matter in the shell (and, 
consequently, of the CR) depends upon the temperature of 

— 3/4 

ambient ISM To, so that vd scales as vd oc T . In Figure 
§we plot v D , defined by @ for T = 10 4 K. 

One can see from comparison of Figures and ^, that 
condition ( fl2| ) is always fulfilled for To = 10 K. It might 
be violated if T ~ 10 6 - 10 7 K, that is if the SB expands 
into the predominantly hot ISM component, but near the 
top the condition ( fL^ ) is always fulfilled even for these high 
temperatures, because Vn oc sin 2 8. 

In the approximation given by (^) the spatial CR den- 
sity perturbations are small and the first term in equation 
( pi] ) is negligible. This means that the drift velocity is al- 
most equal to the velocity of matter vn but opposite to it 
in the direction, so that CR slide through the matter to 
maintain constant spatial density and the time variations of 
their density are only due to the shell expansion. Then the 
equation of continuity reduces to 



(41) d_ 
86 



ncRVii 



A dn C R \ 
1 J 



0. 



ncR 



(44) 



This equation can be integrated with initial condition 
dncR/d9 — at 9 — 0, where Vu = 0, to give 



09 



: ncRVi 



\ucrRJ 



(45) 
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Then the ratio of CR pressure term in the equation ( |10[ ) 
to the gravitational acceleration along the shell is 

Vpgfl _ dpcR dn C R 1 
p 3 g\\ dncR 89 Rgsin9p s 

^P£JL(^3—) 2 S i n 0„io-* sin^ (46) 
p s gR \sm9v D / 

for To = 10 4 K, or maybe ~ 0.1 sintf at To ~ 10 7 K, so that 
the effect the CR pressure gradient has on the sliding of the 
matter from the shell's top is negligible for all interesting 
cases, and, actually, not only at the shell's top, but also for 
9 ~ 1. 

Obviously, the assumption that the flux tube has con- 
stant cross section and the cylindrical geometry implied here 
are not important, because the uniformity of the CR pres- 
sure is achieved due to the very high characteristic drift ve- 
locity vd of the CR compared to the sliding velocity, which 
perturbs the uniformity, and not due to any peculiarities of 
the flux tube geometry. 



9 SUMMARY 

In this paper we consider the effect of the downsliding of the 
matter which takes place in the expanding superbubbles for 
application to the expulsion of the magnetic flux from the 
Galaxy. This expulsion is an important ingredient of the a-fl 
dynamo theory. However, it is shown that even the inclusion 
of the sliding into the calculations of the kinematics of the 
superbubbles, does not enable matter and frozen in flux to 
leave the Galaxy in SBs. One must note, that the impossibil- 
ity of the flux escape from the Galaxy weakens significantly 
the Parker's (1971) argument against the primordial mag- 
netic field, because there is actually no mechanism to expel 
it. 

Some authors (Korpi et al 1999a, Korpi et al 1999b) 
have considered dynamics of the superbubbles in the grav- 
itational field of the Galaxy in greater detail and they also 
find that the matter does not reach terminal velocities larger 
than the escape velocity from the Galaxy. In their simula- 
tions they do observe the development of the SB and its 
blowout from the disk, but at the height of several kilopar- 
secs the velocity of the matter is too small for the matter 
and the field to leave the Galaxy, which agrees with our con- 
clusions. However, they do not comment on its relation to 
the dynamo. 

Other authors (e.g. Hanasz & Lesch 1998, Moss et al 
1999), more directly concerned with dynamos, do not inves- 
tigate the dynamics of escape but merely assume that once 
the magnetic field lines reach the boundary of the disk, they 
are advected away by some mechanism leading to a vacuum 
boundary conditions. They do mention magnetic buoyancy 
as a possible escape mechanism but buoyancy is unlikely to 
be important outside the disk and clearly plays no role in the 
SB escape mechanism, especially when the magnetic fields 
are first amplified from weak seed fields. 

In this paper an analytical formalism for the consider- 
ation of sliding was built on reasonable assumptions, which 
enabled to include the back-reaction of lowering the den- 
sity of the matter on the top of the SB, on the expansion 
of the shell in the radial direction. All the SB, depending 
on their luminosity and conditions in the surrounding ISM 



were demonstrated to fall into two classes: low luminosity 
SBs which are stopped by the Galactic gravitational field 
and fall back, and powerful SBs, which are possibly able to 
reach low density regions at high altitude. Even without slid- 
ing these powerful SBs might be able to expel matter from 
the Galaxy. However, before this happens the shells of these 
high luminosity SBs fragment into separate blobs due to the 
Rayleigh- Taylor instability, which develop when the shells 
start accelerating at the high altitudes, where the density of 
the matter and the outer pressure are very small. These sep- 
arate fragments of the shells continue to move ballistically 
with velocities too small to leave the Galaxy. 

It was shown that our conclusion about the impossi- 
bility of flux to escape in SB explosions does not depend 
essentially upon the details of the inner pressure behavior 
in the SB or the shape of the SB top. The inclusion of the 
CR pressure gradient in the shell also does not influence the 
downsliding because of the high uniformity of CR density in 
the shell caused by the very large diffusion velocity of CR 
in the shell. 

It should be remarked that throughout the paper we 
have assumed that magnetic field lines are tied to the gas. 
If the ionization is the shell is low, then ambipolar diffusion 
(under the influence of the CR pressure gradient perpendic- 
ular to the shell) could allow some small slippage of the field 
lines through the neutral component and, thus, some small 
amount of escape. We do not discuss this possibility in this 
paper. 

All this lead us to the conclusion that if there does exist 
some mechanism responsible for the flux expulsions from the 
Galaxy needed by the a-Sl dynamo this mechanism is not 
superbubbles. The question of whether such a mechanism 
exist remain an open one. 
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